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INTRODUCTION 

In the last decade substantial attention has been 
given to transport of electrons through individual 
atomic point contacts or molecules, commonly referred 
to as quantum junctions (QJ). New experimental re- 
sults for I{V) curves of QJ are promptly followed by 
their numerical modeling which, however, frequently 
lead to results that differ by order of magnitude from 
the measured data pQ. These discrepancies triggered 
interest in the reliability of the often quietly assumed 
approximations. Prime suspects are incorrect atomic 
geometries, inappropriate description of the exchange- 
correlation effects based on ground state methods and 
possible necessity to use a full time-dependent simula- 
tion to arrive to a correct current-carrying stead-state. 

In this short paper we first give a very simple 
derivation of the Landauer formula for a 2-point con- 
ductance of QJ G 2P , based on the uncertainty princi- 
ple. The aim of this is to introduce this central equa- 
tion of quantum transport to a general audience. Next 
we analyse the dynamics of setting up a steady-state 
current in a simple many-electron system and use these 
observations to present physical basis and formal result 
for the 4-point conductance G 4P , rigorously related to 
the non-local conductivity of an extended system con- 
sisting of electrodes and their junction. 

LANDAUER FORMULA AND THE UNCER- 
TAINTY PRINCIPLE 

Let us consider a general QJ. Single-particle quan- 
tum states - wavefunctions that can be occupied with 
electrons and carry current - extend with nonzero 
amplitude from the left electrode through the junc- 
tion into the right electrode. These states, as in every 
metallic system, form a continuum of states of energy 
E e {E min ,E max ), where E mm and E max are the bot- 
tom and the top of the conduction band of the elec- 
trodes in equilibrium. In equilibrium the Fermi energy 
Ep is located somewhere within this interval. Apply- 
ing bias voltage AV between the two electrodes means 
that those states within the continuum that carry cur- 
rent to the right (right-going scattering states) will be 
occupied up to energy fij,, that is, eAV higher than 
the states that carry current to the left (left-going 
scattering states), occupied up to the energy i.e. 
eAV — fiR — Since degenerate energy levels with 
both right- and left- going states occupied carry zero 
total current, the only contribution to the transport 
of charge originates from the energy interval (/zl,/z#), 
occupied only with the right-going states. 

At this point we depart from the traditional deriva- 
tion |2] and make use of the uncertainty principle. 
Right-going electrons occupying states from the energy 
interval ((j>l> ^r) can be put into wavepackets with en- 
ergy uncertainty AE = /i£ — fiR = eAV. The uncer- 



tainty principle then states that there is an uncertainty 
in time At within which we can observe one electron 
passing through the junction 



AEAt 



(1) 



However, there might be only some probability 
T(Ep) < 1, depending in general on the energy 1 , that 
an electron will pass through the junction, out of N 
electrons approaching the junction within time interval 
NAt only fraction NT(Ep) will actually get through. 
Using these observation and the fact that current is 
the number of electrons per time, we find 



I = 2e 
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where we have added a factor 2 to account for the spin 
degeneracy and sum over possibly several degenerate 
right-going states i. This is the celebrated 2 point Lan- 
dauer formula [2] relating the current and the differ- 
ence between electrochemical potentials in electrodes 
Ay, the basic equation of modern mesoscopic physics. 
The quantum of conductance G° = H- = (12.9fcil) — 1 
is obtained for "open" QJ with Ti(Ep) — 1 for only 
one i. Landauer equation represents an important re- 
sult - it relates quantum-mechanical properties of QJ 

- the transmission probability - with the macroscopi- 
cally measured quantity - the conductance. We would 
like to note that the weak "proportional to" sign can 
be made into strong "equal" introducing occupation- 
adapted orthogonal wavepacket basis set. Including the 
mean-field, local-neutrality arguments leads to the 4- 
point conductance, e.g. in ID G 4P = G 2P /(1 - T{E F )) 
relating the current to the induced electrostatic drop in 
potential in the vicinity of QJ, AV 1 [2]- 

DYNAMICS OF ID QUANTUM GAS OF 
ELECTRONS 

While previous derivation gives some hint of the 
time-dynamics in quantum transport, namely that 
based on wavepackets of electrons having significant 
amplitude in a given space for certain time At, the ap- 
proach is really describing a steady state. So we ask 

- how do we make a transition from equilibrium to 
steady non-equilibrium situation in QJ? 

Consider a ID electron gas of density n, fixed by 
£ . At time t = we apply 



Irk 



the Fermi energy Ep 
a localized electric field of the form 



2m 



E' : 



AV 



(6(x + a/2) - 6(x - a/2)), (3) 



where 9() is the unit step function, a the distance on 
which the field is nonzero, modeling the width of QJ, 

1 Assuming smooth dependence of the transmission proba- 
bility on the energy and that the applied bias is small we set 
T(FA Sd TfKd for E <= (n. . u.r V 



1 



and Ay corresponds to the applied voltage. The re- 
sponse of the density and the current can be found 
using the linear response theory In general, the 
non-local conductivity a gives a causal linear relation- 
ship between the total electric field and the current 
density 

J(r, t) = f dt' [ d 3 r'3(f, ¥; t-t')- t') (4) 



For our simple system the nonlocal conductivity, cal- 
culated directly from the continuum of occupied quan- 
tum states, is known analytically 0. Using this within 
the general formula (4) and performing a simple nu- 
merical calculation we find that the current is indeed 
settling to a steady value / = =|-AV with a well de- 
fined relaxation time r (see inset in Fig.l). It is in- 
structive to analyze the dependence of the latter on 
the width of the junction a, while keeping the overall 
bias voltage AV constant. The resulting dependence 
for a gas with Ep = 0.0 7Ha = 2eV (corresponding to 
a gold nanowire) is shown in the Fig.l. The behavior 
of r has a nice physical interpretation: for junctions 




20 40 60 80 100 120 140 160 
a [a.u.] 



Fig. 1. Dependence of the relaxation time r (defined by 
the first local maximum in j(t), see the inset) on the width 
of the junction a. The arrows show the limiting behaviour. 

smaller than the Fermi wavelength Xp, the relaxation 
time is constant and given by the timescale dictated 
by the Fermi energy, for junctions larger than Fermi 
wavelength the relaxation is related to a time it takes 
for an electron with Fermi speed to pass the region of 
the junction. These limiting results are obtainable by 
a detailed analysis of the analytical formulas as well. 

These results demonstrate the mechanism how the 
current is established for non-interacting electrodes, 
which is pertinent to all present ab-initio calculations. 
However, as soon as we include even Hartree inter- 
action, the process changes dramatically. Namely, the 
observed current in Fig.l leads to charge imbalance - 
there are electrons piling up on the right of the junc- 
tion and electrons being depleted on the left of the 
junction. Physically this is almost obvious as any lo- 
calised electric field should be screened out by elec- 
trons. The screening can be overcome only in the limit 
a — > oo, E e — > 0, AV — const when the field becomes 
homogeneous. However, as the Fig. 1 clearly shows, in 
this case the relaxation time becomes infinitely as well! 

The resolution is following. To set a steady current 



on a external homogeneous field of finite magnitude for 
a finite time t < t e . This is to be compared with the 
infinitesimal magnitude and for all positive time in the 
previous paragraph. In the semi-infinite electrode this 
will create an uniform current even for t > t e and there 
is no problem with charge accumulation. Only in the 
region of QJ the charge will pile up on the left and 
deplete on the right. This will lead to a induced field 
E\ localised at the QJ and characterised with some 
induced drop in potential AV 1 . This induced field will 
self-consistently develop to such a form that the cur- 
rent /, established throughout the semi-infinite elec- 
trodes, will be able to pass the junction. 

GENERAL EXPRESSION FOR CONDUC- 
TANCE 

The above ideas can be nicely formalised within the 
framework of the linear response theory. The current 
is given by 

I(x,t)= [ dt'dx'o-(x,x';t-t') (E e (x',t') + E l (x',t')) . 
Jo 

(5) 

Dividing the conductivity into the electrode's, transla- 
tionally invariant part and the part primarily due to 
the junction, <r(x, x') — a H (x — x') + a J (x, x') respec- 
tively, and performing the long-time analysis we arrive 
at the final expression for the steady-state current @] 



I = ^-G 2P AV l = G 4P AV* 



(6) 



where 



a = a H (q = 0,t ^ oo), (7) 
P = -J dq<J J (q,q' =0;t->Oo), (8) 
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where q, q' are wavenumbers, reciprocal to x, x' re- 
spectively. The equation JBJ is expressing the 4-point 
conductance in terms of the limiting character of the 
non-local microscopic conductivity. The utility of this 
formulation lies in a formally straightforward evalua- 
tion of the latter using various ab-initio and, unlike the 
Eq. (J2J, also correlated many-body methods. 
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